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Abstract 

The method of the top quark mass measurement from Mu distri¬ 
bution in t —> hlu decays can be considerably improved if applying 
the technique of moments and proceeding to moments of high degree. 

At the LHC this allows one to reduce the systematic and statistical 
errors of the top mass by more than a factor of two. 

The precision measurement of the top quark mass is necessary for testing 
of the SM and/or selecting the probable scenario of its extension. The current 
accuracy of the top mass, attained at Tevatron, is A M t = 2.9 GeV, and 
a further progress is expected mainly at LHC and a future International 
linear collider. At the LHC the accuracy of measurement of the top mass is 
anticipated at the level of 1- 2 GeV and in view of large statistics the main 
difficulties are expected from systematic errors. An analysis of Tj shows 
that the most promising method from the point of view of optimization of 
the errors is based on an investigation of the invariant-mass distribution of 
Jfi/j + f system produced in decays t —> bW —> Mu with 6-jet formation and 
the J/0 being a member of the 6-jet. The evaluation made by Monte-Carlo 
modelling gives approximately 0.7 GeV for systematic error of the top mass 
measured by this method and of about 1 GeV for statistical error for 4 years 
of LHC operation j2]. In this paper we discuss an opportunity to improve 
this result by means of optimization of the data handling within the same 
experiment. 

*Talk given at 12th Lomonosov Conference on Elementary Particle Physics, Moscow, 
Russia, August 25-31, 2005. 
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So, following we consider the processes 

qq, gg —> tt —> bW bW —> 6 f?z/ 691^2 =>• { 6 -jet + £} + {3 jets} , (1) 

with the 6 -jet and the isolated lepton £ = {e, yu} coming from one t quark, 
and the remaining three jets coming from another t quark. In the experiment 
the above states are registered and the invariant-mass distribution of J /0 + £ 
system is measured. The Mj/^n distributions can then be compared with a 
template of shapes parametrized by the top mass and thus M t can be fitted 
| 5 - One could note, however, that the invariant-mass distribution of the b+£ 
system inevitably emerges at a certain stage of this analysis. Therefore 011 
the equal rights one can carry out the fitting in terms of the data converted 
to the form of Mu distribution. In what follows we adhere to the latter 
option and hence consider that (as if) a distribution F(q) = a~ l da/dq is 
measured, where a is the cross-section of the process {TJ) and q = M m is the 
reconstructed invariant mass of the 6 + £ system. 

By proceeding in this manner, we simulate the data under the follow¬ 
ing suppositions. First we suppose that there is a satisfactory method for 
extracting signal from the data. (Actually this means the existence of a sat¬ 
isfactory model for background processes that survive after setting of kine¬ 
matic cuts [21.) Secondly we describe the signal in the Born approximation, 
identifying the 6 -jet with the 6 quark. Finally, on the basis of the results of 
|5, we disregard the effects of the finite width of the top quark. The lat¬ 
ter assumption means that o~ x d<jjdq is equal to F^} dr y M / dq, where Fuv 
is a partial width of the decay t —> b£v. (Thus the distribution F becomes 
process-independent.) 

At Fig.^a) the distribution F(q) is represented for various values of the 
top mass. We see from the figure that in the absolute value the distribution 
is most sensitive to M t in a region located between the maximum of the dis¬ 
tribution and a large-g tail where F(q) is almost vanishing. So in this region 
of g-variable one could expect the highest accuracy for the M t determina¬ 
tion. A practical way to enhance the role of this region is to proceed to high 
moments over the distribution, 



( 2 ) 


Here M is a technical parameter fixed close to M tl and the normalization 
of F(q) is adjusted so as to satisfy the equality (1) = 1. From Fig.[T{b) it is 
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obvious that with increasing n the distribution q n F(q) is concentrating more 
and more in the above intermediate region, increasing thus the sensitivity of 
(q n ) relative to M t . 

So, by basing on the above observation, we define the experimentally 
measured top-quark mass as a solution to the equation 


<?">«P = ( q n ) . 


( 3 ) 


Here the moment in the l.li.s. is determined on the basis of the experimental 
data, and that in the r.h.s. on the basis of the theoretical distribution, which 
depends on the parameter M t . Let, for a given n, a solution to © be 
M t = Mt(n)- Then the error of the solution is 

^ 

Our aim is to estimate A M t i n ) and hnd an optimal value of n which would 
minimize A M t r n y Since in view of (J2J) the derivative d(q n )/dM t is known, the 
problem is reduced to the estimation of the error A (q n ) exp or its components, 
the statistical and systematic experimental errors. 

The statistical error can be estimated, in almost model-independent way, 
by means of the formula |3, 



A stat (g n ) eX p 



( 5 ) 


Here N is the representative sample of events, and ( q 2n ) is a theoretically 
calculated moment. For the consistent determination of the systematic error 
the application of a proper MC event generator is needed. The analysis [2j 
made with the aid of the PYTHIA and HERVIG events generators showed 
that the major uncertainties in the Mj/^ e distribution were caused by (i) the 
uncertainties in the 6-quark fragmentation and (ii) the background processes. 
It is obvious that at solving the inverse problem, the reconstructing the My 
distribution from the one, the errors should be of the same origin. 

With this in mind we consider first the error resulting from the uncertainty 
in the 6-quark fragmentation. For brevity we call this the type I error. At 
the level of the My distribution it appears as an uncertainty A q in the q 
variable. Let us suppose that A q is sufficiently small, and let us neglect 
nonlinear effects. Then we have 

fM 

A sys V)ex P = J o dq [q n F(q)]' A q . (6) 
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Here the prime means the derivative with respect to q. The A q is determined 
by the relation A q = rq with r is a coefficient, r ~ ^AE b /E h} where E b is the 
energy of the b quark in the laboratory frame and A E b is its error. In fact 
the mentioned relation for A q follows from the observation that q 2 actually 
is a doubled scalar product of 4-momenta of the b quark and of the lepton l 
whose momentum is considered precisely determined. 

The systematic error arising after subtraction of the background processes— 
we call it the type II error—appears in the absolute value of the distribution 
function. Denote the corresponding error by SF, then we have 

r M 

A sysI V) exp = j o d q q n SF(q). (7) 

The function SF(q) should be vanishing at the boundaries of the phase space. 
We also assume that SF only once changes the sign when q runs the values. 
The simplest function satisfying to these requirements is a polynomial, SF = 
hq(q—M/2)(q—M) with parameter h describing the amplitude of the error. 

Now we turn to numerical results. We assign the following values for 
the parameters with global meaning: M w = 80.4 GeV, Y w = 2.1 GeV, 
M t = M = 175 GeV. The parameters N, r, h take the values depending on the 
conditions of observation. Since in |2J the Mm distribution was determined 
at N = 4000 (with kinematic cuts and for 4 years of LHC operation), we set 
this value for N , as well. The parameters r and h may be fixed by basing 
on the properties of the Mj^ i distribution. We omit here the respective 
discussion and show only the result: r ~ 0.005, h ~ 1.7 x 10~ 10 GeV -4 
J2|. Note that the same estimation for r follows from the above mentioned 
formula r ~ -A E b /E b when taking into consideration the 1%-precision of 
the determination of the energy of the b jets expected at LHC [TJ. 

Now, as we know the parameters, we can calculate the errors of the 
top mass. At Fig. 12(a) we show the behavior of A stat M i(n ) versus degree 
n of the moments. Fig. 0b) shows the behavior of the summed-quadrature 
systematic error A sys M t ( n \. The dashed lines in both figures represent the 
errors obtained with usage of effective moments determined by way of in¬ 
troducing a cut-off A„ instead of M in formulas 0 , @ and 0 ; the aim 
of introducing the cut-off is to isolate the tail of q n F(q) with an unphysi¬ 
cal peak arising at large n, see Fig. 0b). Having the errors A stat M t ( n ) and 
A sys ilfy( n ), it is easily shown that near n — 15 the total error reaches the 
minimum AM t ( 15 ) = 0.48(0.42) GeV, with A stat M t ( 15 ) = 0.41(0.40) GeV 
and A sys M t ( 15 ) = 0.24(0.08) GeV. Here we put into the brackets the errors 
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Figure 1: The distribution F(q)=T 1 dT/dg, q = Mm, at M t = 170, 175, 180 GeV (a), 
and the shape of q n F(q) at n = 1, 5, 15, 40 (M t = 175 GeV), arbitrary normalization (b). 



calculated by the effective-moment method. It can be shown by the naive 
counting method |3J that the theoretical errors will not spoil the above accu¬ 
racy provided that the M b( i distribution will be calculated within the one-loop 
electroweak and two-loop QCD precision. 

In summary, we discussed the potentialities inherent in the technique of 
moments and demonstrated that at the LHC the proceeding to moments of 
high degree allows one to reduce the total error of M t to the level of about 
500 MeV. 
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